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ABSTRACT 



The GOES M2-class solar flare, SOL2010-06-12T00:57, was modest in many 
respects yet exhibited remarkable acceleration of energetic particles. The flare 
produced an ~50 s impulsive burst of hard X- and 7-ray emission up to at least 
400 MeV observed by the Fermi GBM and LAT experiments. The remarkably 
similar hard X-ray and high-energy 7-ray time profiles suggest that most of the 
particles were accelerated to energies ^300 MeV with a delay of ~10 s from mildly 
relativistic electrons, but some reached these energies in as little as ~3 s. The 
7-ray line fiuence from this flare was about ten times higher than that typically 
observed from this modest GOES class of X-ray flare. There is no evidence 
for time-extended >100 MeV emission as has been found for other flares with 
high-energy 7 rays. 

Subject headings: Acceleration of particles — Sun: flares — Sun: particle emis- 
sion — Sun: X-rays, gamma rays 



1. Introduction 



The Sun is capable of accelerating electrons and ions to relativistic energies on time 
scales as short as a few seconds during solar flares. This conclusion has been reached based 
on observations of the X rays, microwaves, 7 rays, and neutrons produced wh e n the flare- 



accel erated particles interact in the solar atmosphere (IForrest fc ChuppI Il983l : iKane et al. 



19861 ). The first reported observation of 7 rays with energies above 10 MeV w as made with 



the S olar Maximum Mission (SMM) spectrometer during the 1981 June 21 flare (jChupp et al. 



19821 ). Most of the 7-ray emission occurred within an ~70 s period and was followed minutes 



later by detection of high-energy neutrons at the spacecraft. Although it was clear from the 
neutron timing observations that protons were accelerated to energies in excess of 100 MeV, 
it was not possible from 7-ray spectroscopic studies to conclude that the protons reached the 
energies a bove ~300 MeV nece ssary to produce the characteristic spectrum from pion-decay 
radiation (IMurphy et al.l 119871 ). 
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Forrest et al.l (119851 Il986l ) provided clear spectroscopic evidence for pion-decay emission 
during both the prompt and delayed emission phases of the 1982 June 3 flare. This emission 
in the second phase was confirmed by the presence of the 0.511-MeV annihilation line from 
the decay of positi vely charged pions that had a time profile similar to the high-energy 7-rays 



(IShare et al. 1119831). 



1 



Since the era of those early measurements, improved spectrometers have detected the 
presence of pion-decay emissio n in several more flares and in two cases have observed 7-ray 



emission up to at least 1 GeV ( Kanbach et al. 19931: 



Kuznetsov eraPboilh . iRvanI J2000h and IChupp k. Ryad J2009h have reviewed these and 



Akimov et al. 1996: Vilmer et al. 2003 



other examples of high-energy 7-ray emission in flares. Most of these high-energy emissions 
have been observed over tens of minutes to hours leading to the designation of a class known 
as 'Long Duration Gamma-Ray Flares'. Several processes have b een suggested to explain 
such particle acceleration to high energies in the solar environment (IE . 



Petrosian et al. 1994: Park et al. 1997 



Rvanlboool : [AschwandenlboO? 



lison fc Ramatv 



Chupp fc Rvanll2009h . 



1985 



Typically such high-energy 7-ray emission has been associated with intense soft X-ray 
flares recorded by (JOE'S" having peak powers exceeding 10~^ W m~^ (X class). It is of interest 
to determine whether less energetic solar flares have the capability of accelerating electrons 
to energies ^100 MeV and protons to energies >300 MeV. Out of 24 flares observed at > 10 
MeV energies by the 7-ray spectrometer on SMM only 3 had GOES classifications of M5 or 
less ( IVestrand et al.lll999l ). Similarly, most flares emitting 7-ray lines have been associated 
with intense soft X-ray emission: out of 65 7-ray line flares observed by SMM only three have 
been associated with flares having a GOES classification of M5 or less. Before detection of 
this flare, the Ramaty High Energy Solar Spect roscopic Imager ( RHESST) observed only one 
nuclear-line flare out of 20 below the M5 class ( IShih et al.ll2009l ) and that was characterized 
as C9.7. The smallest GOES class flare for which detection of nuclear 7-rays has been 
claimed was the C7 flare observed by the COMPTEL Compton telescope on the Gompton 
Gamma-Ray Observatory (GGRO) (lYoung et al.l 1200 ll ). 



With the launch of the Fermi mission in 2008, it is now possible to make the high- 
sensitivity measurements necessary to detect > 30 MeV 7 rays in the weakes t flares. Fermi 



is co mprised of two instruments: the Gamma-ray Burst Monitor (GBM) (IMeegan et al. 
2OO9I ) sensitive from ~ 8 keV up to 40 MeV and co vering the energy ba nd of nuclear 7-ray 



line emission; and the Large Area Telescope (LAT) ( lAtwood et al.l 120091 ) operating from 20 
MeV to more than 300 GeV and covering most of the pion-decay emission energy range. In 
its normal sky survey mode the LAT observes the Sun for only ~ 35 min every 3 hours. It 
was therefore fortunate to observe the M2 GOES-class flare on 2010 June 12 (SOL2010-06- 
12T00:57). This flare produced nuclear line emission about an order of magnitude higher 
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than typical for this magnitude soft X-ray flare. Although the flare lasted only about 2 
min, it appears to have accelerated electrons and/or protons to energies >300 MeV based 
on detection of 7 rays with energies >100 MeV. One of the key features of this observation 
is that the high-energy emission is delayed by ~ 5 — 10 s relative to the 500 - 1000 keV 
bremsstrahlung. This suggests that the acceleration of hundreds of keV electrons and hun- 
dreds of MeV electrons and/or protons in the chromosphere took place within 10 s of each 
other. A weak solar energetic particle (SEP) event observed by GOiJS' followed this flare. 

In the next two sections we discuss the capabilities of the GBM and LAT for detecting 
7-rays from solar flares. These are followed by §4 describing the observations including 
spectroscopic studies in the nuclear energy range with the GBM and in the high-energy 
domain with the LAT. As most high-energy 7-ray flares previously observed exhibited an 
extended emission phase, we then discuss a search for such emission in this flare in the next 
section. In the final section, §5, we summarize our conclusions and discuss the results. 



2. GBM Capability for Flare Observations 

The Gamma-ray Burst Monitor (GBM) was designed to observe gamma-ray bursts but 
has useful capabilities for other sources such as SGRs (soft 7-ray repeaters), pulsars. X-ray 
binaries usin g the Earth occultat ion technique. Terrestrial Gamma-ray Flashes (TGFs), and 



solar flares (IMeegan et al.ll2009l ). The GBM is comprised of twelve sodium iodide (Nal) 
detectors measuring the energy range from ~8 keV to 1 MeV, and two bismuth germanate 
(BGO) detectors covering the range from ~200 keV to 40 MeV. The detectors are arranged 
to collectively view the entire sky not blocked by the Earth. 

In this paper we use the BGO detector viewing the Sun to observe the nuclear line 
and continuum emission from solar flare SOL2010-06-12T00:57. Each BGO detector is a 
cylinder of length and diameter 12.7 cm, viewed by a photomultiplier tube at each end. Its 
effective area for detecting photons ranges from 160 cm^ to 200 cm^, depending on energy 
and direction of incidence. Importantly for measuring solar flare nuclear lines, the high Z of 
bismuth, the high density of BGO, and the large volume of the GBM BGO detectors result 
in a high probability for absorbing the full incident photon energy: ~ 67% at 1 MeV, 50% 
at 3 MeV and 40% at 10 MeV. 

Several types of data are produced by the GBM. There are two temporally-binned types. 
Continuous Time (CTIME) and Continuous Spectroscopy (CSPEC), as well as Time- Tagged 
Events (TTE). For intense solar flares, TTE data are likely to be lost due to a bandwidth 
limit between the GBM and the spacecraft solid-state recorders. Therefore for solar-flare 
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spectral analysis the appropriate data type is CSPEC. CSPEC has 128 quasi- logarithmic 
energy bins and 4.096 s temporal resolution, which is improved to 1.024 s when the flight 
software triggers on a statistically significant rate increase. 

The detector performance was measured before launch using X-rays from the BESSY 
synchrotron radiation facility, radioactive isotopes including the ^^C line at 4.43 MeV (from 
an ^^^Am/^Be source), and 7-rays up to 17.6 MeV produced via (p,7) reactions using a 
Van de Graaff accelerator. Th ese data have been used to calibrate th e gain, resolution. 



and response of the detectors (IBissaldi et al.l l2009t iMeegan et al.l l2009l ) . The automatic 



gain software measures and adjusts the gains of the BGO detectors using the 2.223 MeV 
line present in the background from capture of cosmic-ray produced neutrons on H in the 
spacecraft's propellant. In order to improve the gain solution for analysis of the June 12 
flare data we measured the centroid of the solar neutron-capture line. This resulted in a 
1% gain change that is small relative to the ~ 8% resolution of the GBM at the 2.223 MeV 
line. The best-fit energy of the positron annihilation line (0.511 MeV) from the flare was 
found to be 0.530 ± 0.007 MeV, in disagreement with fits of the instrumental background line. 
However, both lines are superimposed on strong continua making it difficult to determine the 
lines' peak position. Further investigation of the BGO calibration below 1 Me V is planned. 



parti cularly since TGF spectra also show a shifted positron annihilation line (iBriggs et al. 



20111) 



Spectral fitting is performed using the forward-folding technique with an assumed parametrized 
photon spectrum folded through a detector response matrix (DRM) to produce a counts 
spectrum; the parameters are adjusted to obtain the best fit to the observed counts spec- 
trum. The DRM is based on Geant4 simulations of the GBM detectors and mass model of 



the s atellite and has been validated by comparison with the calibration data (IHoover et al. 



20081). 



3. LAT Capability for Flare Observations 



The Large Area Telescope (LAT) is a pair-c onversion telescope designed to detect 
gamma rays from 20 MeV up to more than 300 GeV (lAtwood et al.ll2009l ). It is made up of a 
4x4 matrix of identical towers, each one comprised of a tracker with layers of Silicon Strip 
Detectors (SSD) alternating with foils of high-Z converter (tungsten), and a calorimeter with 
logs of Csl arranged in a 'hodoscopic' configuration so that the energy deposition is imaged 
in three dimensions. The array of towers is surrounded by an Anti-Coincidence Detector 
(ACD) made up of 89 plastic scintillator tiles with a 5 x 5 array on the top and 16 tiles 
on each of the 4 sides. To reduce the impact of the self-veto events caused by calorimeter 
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back-splash, each ACD tile is typically <1000 cm^, depending on its position in the array. 
A 7-ray passes through the ACD with small probability of interaction and can convert into 
an electron-positron pair that is tracked in the SSD. The energy of a photon below about 
100 MeV can be estimated by multiple-Coulomb scattering of the electrons in the tungsten 
converters and SSD. Energies of higher-energy photons are measured by 'total' absorption 
or modeling the shower profile in the calorimeter. In the standard Fermi sky-survey mode 
the spacecraft rocks 50° north and south in celestial declination from the zenith so that each 
region of the sky is viewed for ~30-40 min every two orbits; therefore the LAT's large aper- 
ture (2.4 sr) and effective area provides the capability to sensitively monitor solar activity 
with a duty cycle of 15-20%. 

The ACD is used to reject the large background of charged cosmic-radiation and secon- 
daries from the spacecraft and Earth's atmosphere. The threshold for this veto is nominally 
set at 45% of the amplitude of a minimizing ionizing singly-charged particle traversing a tile 
in the ACD (i.e. set to ~800 keV). If this threshold is exceeded and if the ACD tile hit 
was adjacent to a tracker tower that caused the event to trigger, the event will be vetoed 
unless a sufficiently high energy is deposited in the calorimeter (>100 MeV in one or more 
crystals). This veto inhibition ensures that backsplash does not cause very high-energy 7 
rays to self-veto. Of the vetoed events, 2% are telemetered to the ground for diagnostic 
purposes. 

Solar flares can emit intense fluxes of tens of kcV X rays. About 20% of the X rays at 
these energies can penetrate the thermal blanket and micrometeoroid shield and reach the 
ACD, depositing part of their energy in the illuminated tile. 

Several of these photons can arrive within the 0.4 fis anti-coincidence veto shaping time 
(pulse pile-up) to yield a high total-energy loss. It is possible in the largest of flares that this 
energy loss could exceed the 800 kcV veto threshold and information about a valid photon 
event would not be transmitted to the ground. None of the solar flares detected to date have 
been sufficiently intense to cause this data loss to occur. 

Pulse pile-up of hard X rays from some flares has, however, affected the ground analysis 
of LAT data, in the classification of events as 7 rays or background. ACD tile hits are 
registered (i.e. included in an event) whenever there is >100 keV energy deposition integrated 
over the ~3 /is peaking time of the front-end amphfier (as opposed to the 800 keV threshold 
for an ACD veto). Pulse pile- up during even modest flares generates such tile hits. In order 
to achieve the highest sensitivity for studying celestial sources, the current LAT ground- 
processing software rejects 7-ray events with a high ratio between the number of ACD tile 
hits (i.e., the number of tiles with energy deposition >100 kcV) and the energy measured by 
the calorimeter. Thus, 7 rays arriving during relatively intense flares have a high likelihood 



- 9 - 



of being rejected from the standard LAT data products. This is the case for the 2010 June 
12 flare as we will discuss in §4.2 below. 



As a consequence of this feature, the LAT Low- Energy (LLE) technique (jPelassa et al. 



20101 ) was adapted for analyzing data during this flare. LLE event selection uses less dis- 
criminating criteria than the standard ground-processing and is not affected by ACD tile hits 
>100 keV. The primary requirement is that the candidate 7-ray event have at least one track 
and a reconstructed energy larger than 30 MeV. Due to these less discriminating criteria the 
off-axis angle for LLE events can be as large as ~80° compared with ~68° in the standard 
LAT data products. In addition, only 7-rays whose estimated arrival directions were within 
20° of the Sun are included for analysis. The 20° angular restriction was chosen based on 
Monte Carlo simulations of the LAT point-spread function at these energies. A time series 
of events is then constructed from which background intervals on either side of the flare are 
deflned and a linear or quadratic interpolation is used to estimate the background during the 
flare. A Detector Response Matrix (DRM) for the solar location during the flare is created 
using a custom Monte Carlo simulation. By passing candidate photon models through the 
DRM, we then flt the background-subtracted data using a version of rmflt 3.4|l |, customized 
for the speciflc solar flare, and the OSPE}(§ | analysis packages. 

Using a Monte Carlo analysis, we estimate that the energy resolution for a source at ~75° 
off axis in our LLE analyses is about 40%. We also estimate ~15% and ~10% systematic 
uncertainties in the LLE fluxes >30 and 100 MeV, respectively, based on studies of the Vela 
pulsar. The uncertainty >100 MeV is consistent with that found using the standard Fermi 



analysis tools applied to the pulsar (lAbdo et al.l 120101 ). 



4. June 12 Flare Observations 

The GOES M2.0 class X-ray flare commenced with some low-level activity on 2010 
June 12 at 00:30 UT. Although both RHESSI and Fermi observed the flare, here we only 
discuss the Fermi observations because RHESSI was not pointed at the Sun during the 
flare. The flare occurred with Fermi in sunlight and during a relatively low-background 
portion of its orbit. As seen in Figure [1^ the 10 — 26 keV emission recorded by the GBM 
Nal detectors commenced around 00:55 UT and rose precipitously about 40 s later; for 
comparison we also plot the GOES 0.5 - 4 A proflle and note that this emission is dominated 



iR.S. Mallozzi, R.D. Preece, & M.S. Briggs, "RMFIT, A Lightcurve and Spectral Analysis Tool", Robert 
D. Preece, University of Alabama in Huntsville, (2008): |http: / /fermi.gsfc. nasa.gov/ssc/data/analysis/user /] 

^SolarSoft: http:/ /www.lmsal.com/solarsoft/ 
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by 3 keV thermal photons as is reflected in its slower rise and extended tail. The GBM burst 
algorithm triggered on the flare at 00:55:05.64 UT and put the instrument in a high-data 
rate mode for the next 10 min. The 100 — 300 keV time proflle observed by the GBM's solar 
facing Nal detector is also plotted in Figure [T^. It is clear that the emission peaks more 
sharply and ends sooner at higher X-ray energies. Most of the emission observed above 100 
keV occurred within an ~50 s interval. The event as viewed by GOES 1 - 8 A channel 
ended about 01:30 UT. The flare originated from active region (AR) 11081 at approximately 
N23°W43°. White light emission observed b y the Helioseismic and Magnet ic Imager (HMI) 



on the Solar Dynamics Observatory (SDO) ( iMartmez Oliveros et al.ll201ll ) in a single 45 s 



exposure, consistent in time with the hard X-ray emission, reveals two compact footpoints 
about 10^ km apart. 

There was no evidence for signiflcant dead time and/or pulse pile- up effects in the 
GBM BGO detector facing in the solar direction. Photons with energies up to ~8 MeV 
were detected by the GBM during the 50 s peak. At the time of the flare the spacecraft was 
rocking 50° to the south so that the Sun was 76° off axis, close to the edge of the fleld-of-view 
(FOV) for LLE studies, and the Earth's horizon was entering the FOV. The accompanying 
hard X-ray emission from the flare was detected in the LAT's ACD. In Figure [T]d we plot the 
average number of ACD tile hits as a function of time. As discussed in the previous section, 
pulse pile-up from tens of keV hard X-rays exceeded the 100 keV ACD hit threshold. This is 
reflected in the broad peak with a maximum near 00:57 UT that has a shape similar to the 
11-26 keV emission observed by the GBM Nal detector. The impulsive peak in the ACD rate 
is also similar to that observed between 100 and 300 keV by the GBM Nal detector. There 
is no evidence for an increase in the number of ACD vetoed events in orbit during the flare, 
indicating that pulse-pileup from hard X-rays did not exceed the 800 keV veto threshold. 
Thus, the overall valid event rate transmitted to the ground for processing was not affected 
by the ACD response. However, as shown by the red curve in Figure ^ there is no evidence 
for the flare in the well-screened standard LAT data products (show n in the figure are the 



events belonging to the "transient" event class, lAtwood et al.ll2009l ). If anything, we see a 
deficit of events in the standard analysis light curve which is a consequence of the high ACD 
hit rate >100 keV, shown in Figure [Do, that caused a significant increase in the number of 
events that failed the standard quality cut. 

It is important to establish convincingly that any high-energy emission observed by the 
LAT originated at the fiare site and was not due to artifacts from the high rates encountered 
by the instrument. The most compelling evidence for the solar fiare origin is the map of 
events relative to the position of the Sun. Plotted in Figure [2^ is the distribution of LLE- 
selected events with energies >30 MeV accumulated 30 s before and after the fiare. The 
distribution is affected by the instrument field of view and by the removal of events near the 
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Fig. 1. — Time histories related to the 2010 June 12 solar flare, a) GOES 0.5 - 4 A rates, 
and GBM Nal 11 - 26 keV and 100 - 300 keV relative rates; b) LAT ACD hit rate >100 keV 
containing contributions from background, >100 keV solar flare X rays (impulsive peak) and 
pulse pile up from lO's of keV solar X rays following the Nal 11 - 26 keV profile in la); and 
c) LLE and LAT Transient Class event rates. 
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Earth's horizon. We reduce the contamination from the bright Earth hmb selecting only the 
events with a reconstructed zenith angle less than 100°, but, due to the large point spread 
function at low energy, some residual contamination near the horizon may be due to this 
limb brightening. 

In Figure [2}3 we plot the distribution accumulated from 00:55:40 to 00:56:30 UT dur- 
ing the flare. We note that the excess near the Sun's position is biased due to its loca- 
tion near the edge of the LAT's field of view. This is known as the "fisheye" effect, i.e., 
the ten dency of events to be reconstructed with directions closer to the z axis than they 



should ( iThompson et al.lll993[ ). This bias is relevant only for off-axis events, and it is par- 
ticularly evident at low energies. The position of the centroid of the counts is shown by the 
filled star (at R.A.=84°.39, Dec. =18°. 79, J2000). It is important to note, that in the routine 
LAT analyses the contribution to the overall exposure and photon counts from far off-axis 
events is negligible. In addition to this, owing to the scanning mode of observations with 
the LAT, persistent sources will have the great majority of their exposure at much smaller 
inclination angles. 

In Figure [2t we plot the background-subtracted distribution observed during the flare. 
There appears to be significant 7-ray emission consistent with the reconstructed position 
of the Sun up to energies of ~400 MeV. We have confirmed the observed 6° shift in the 
position of the centroid from the solar position using a Monte Carlo study. Using the shift 
and identifying a < 20° acceptance angle (dashed curve in the figure) for LLE events >30 
MeV maximizes the signal/background. Events that meet these criteria make up the LLE 
rate (black curve) plotted in Figure [It. The overall rate decreased, especially after 01:00 
UT because the defined aperture approached the exclusion region near the Earth's limb and 
many events were consequently rejected. The >30 MeV LLE time profile during the flare is 
similar to the 100-300 keV Nal time history. 



4.1. GBM Spectroscopic Studies 

In order to obtain background-subtracted spectra we used GBM/BGO spectra accu- 
mulated just before the flare and after 4-min following the flare for background in order to 
avoid times when there was significant 15 - 50 keV hard X-ray emission. The background- 
subtracted GBM counts spectrum accumulated over a 50-sec period (00:55:40-00:56:30 UT) 
during the impulsive phase of the flare is shown in Figure |3l Although the flare was relatively 
weak and the BGO detector has only moderate energy resolution, line features are clearly 
evident in the spectrum. 
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Fig. 2. — Angular distribution of >30 MeV 7 rays relative to the position of the Sun detected 
using the LLE analysis. The open star shows the location of the Sun which is close to the is 
close to the edge of the LAT field of view. The filled star shows the shifted position due to 
the fisheye effect discussed in the text. The dotted curve depicts the 20° region containing 
the events used in the analysis, a) Average of background distribution taken 30 s before and 
after the flare, b) Angular distribution of events observed between 00:55:40 and 00:56:30 
UT. c) Difference between angular distribution observed during the flare and the average 
background distribution. The distributions have been smoothed (with a Gaussian kernel) to 
reduce statistical fluctuations. 
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Fig. 3. — Background-subtracted BGO count spectrum accumulated between 00:55:40 and 
00:56:30 UT. The spectrum has been fit by a simple power law (blue histogram) and a 
flatter power-law with exponential cut off (light green histogram) to model bremsstrahlung 
by electrons (see text), and with a nuclear de-excitation component (purple) plus 0.511 and 
2.223 MeV lines (dashed and dark green), and a pion-decay (orange) component. 
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The spectrum has been fit with a photon model consisting of several components. We 
find that the electron-bremsstrahlung component has the shape of power law at low energies 
that hardens above several hundred keV and then rolls over in the MeV range. In this paper 
we fit the bremsstrahlung spectrum with the sum of a low-energy power law and a flat power 
law with exponential cutoff in the MeV range (blue and light green curves, respectively, in 
Figure [3l). Such a complex shape has been observed before in the spectra of several flares 
detected by 5'MM and RHESSI (Share & Murphy, in preparation). iMcTiernan fc Petrosian 
( I1990I ) have found that the magnitude of the hardening above a few hundred keV observed in 
some flares is larger than that e xpected for an electron spectrum following a single power law. 
Rieger &: Marschhauserl (jl99ll ) describe flattening in the MeV range followed by roll overs 
above several tens of MeV in the spectra of soin e flares. Such features cannot be explained 
by transpor t effects alone (jPetrosian et al.lll994l ) and must be produced by the acceleration 
mechanism. iPark et al.l (Il997l ) show that these features can be explained by models based on 
stochastic acceleration by turbulence once loss mechanisms are properly included. Whether 
such models can explain the 2010 June 12 bremsstrahlung spectrum requires more study. 

The nuclear de-excitation lines and continua are represented by a template based on a 
detailed study of nuclear gamma-ray pro duction from accelera ted-particle interactions with 
elements found in the solar atmosphere ( iMurphy et al.l 120091 ). Such templates depend on 
the assumed ambient composition and accelerated-particle composition and spectrum. The 
BGO spectral data for this flare are inadequate to distinguish among templates derived for 
different ambient abundances, particle spectra and angular distributions. For this reason 
we used a arbitrary templates based on earlier studies (Share & Murphy, in preparation). 
The ac celerated part icles were assumed to interact in a thick target with a coronal compo- 
sition (jReamed 119951 ) ■ but with '^He/H = 0.1. We do not take into account transport effects 
but instead assume that the accelerated particles have a power-law differential spectrum 
{dN/dE oc E^ with /3 = —4), coronal elemental abundance (but with an accelerated a/p 
ratio of 0.2), and an angular distribution that is isotropic in the downward hemisphere. As 
we mentioned above our results are not sensitive to the assumed compositions. The pho- 
ton model also includes Gaussians representing the 0.511 and 2.223 MeV positron-electron 
annihilation and neutron-capture lines, respectively, and a pion-decay spectral component. 

We present the best-fltting spectral parameters in Table 1, along with estimates of their 
la uncertainties. We list fluences obtained by integrating over the 50-s time period, except 
for the 0.511 and 2.223 MeV lines. The 0.511 MeV line originates from radioactive decays 
( iKozlovsky et al.l Il987l . 120041 ) with half lives that can extend up to hours, in addition to a 
prompt component from positron annihilation following decay of positively charged pions. 
The neutron-capture line is also delayed because of the time required for the neutrons to slow 
down in the solar atmosphere and photosphere and be captured on H, forming deuterium 
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with the release of 2.223 MeV of binding energy from the mass excess. For these two hnes 
we have fit the background-subtracted spectrum integrated over a total of 250 s. 

The low-energy power law (PLl) is well defined by the fit. The higher-energy power law 
(PL2) with exponential rollover is not as well defined because it competes with the nuclear 
de-excitation line spectrum. However our studies indicate that this component is required 
to provide an acceptable fit to the June 12 flare spectrum. With this component present, 
the fit is acceptable at the 10% confidence level based on the statistic; without it, the fit 
is rejected (0.001 % confidence level). There is no evidence for pion-decay emission in the 
BGO spectrum. As can be seen in Table 1, assuming an hadronic origin, the LAT-measured 
pion-decay fiuence that we discuss in the next section is about a factor of ten lower than 
detectable by the GBM BGO detector. 



4.2. Combined LAT and GBM Spectroscopic Studies 

We have obtained a background-corrected LAT count spectrum >30 MeV accumulated 
during the 50-s period 00:55:40-00:56:30 UT using the LLE data plotted in Figure [It. This 
spectrum revealed flare emission up to an energy of ~400 MeV. The fundamental question 
is: what is the origin of this emission? The nuclear line emission observed with the GBM 
implies the presence of accelerated ions up to at least 50 MeV nucleon^^. It is possible that 
the fiare-accelerated proton spectrum extended up to the ~300 MeV threshold for pion pro- 
duction. Alternatively, it is possible that the LAT emission is from electron bremsstrahlung, 
either from an extension to high energies of the electron spectrum producing the X-ray 
bremsstrahlung observed in the GBM or from an additional hard electron component. One 
possible way to resolve this ambiguity is to jointly fit the GBM and LAT spectra assuming 
different origins for the LAT emission. 

In Figure m we plot the background-subtracted photon spectrum from 0.3 to 400 MeV 
including both the GBM and LAT data. We made two fits to the joint data sets, one 
assuming that the observed LAT emission was from pion-decay radiation (top panel) and 
the other assuming that it was from a hard power-law spectrum of electron bremsstrahlung 
(bottom panel). Based on the statistical quality of the fits to the LAT spectrum we cannot 
distinguish between the two emission models. In addition we cannot constrain the origin of 
the emission for this event by extrapolating the models into the GBM energy range; however 
we note that for a stronger flare we might be able to rule out a power-law model. The choice 
of model also has only a small effect on the parameters listed in Table 1 derived from fits to 
the 0.3 to 40 MeV GBM data. Also plotted in the figure are the extensions into the LAT 
energy range of the power-law and cutoff power-law components derived from the fits to the 
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Fig. 4. — Combined GBM/LAT photon spectrum accumulated between 00:55:40 and 
00:56:30 showing the best total fit using the same components as in Figure 3 plus an addi- 
tional component for the LAT emission. The upper panel shows a pion-decay fit to the LAT 
spectrum; alternatively the lower panel shows a power-law fit, presumedly representing a 
third electron bremsstrahlung component. Note that because this is a photon representation 
the lines are plotted at their intrinsic resolution and appear to be more significant than they 
really are. 
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GBM data. The intensities of these components fall at least an order of magnitude below 
the LAT measurements and therefore do not make a significant contribution to the solar 
emission observed by the LAT. 

Even though we cannot statistically distinguish between a pion-decay or electron-bremsstrahlung 
origin for the observed LAT emission, we can obtain the best-fitting parameters for these 
components. If the LAT emission is from electron bremsstrahlung, we have shown that it 
cannot be a simple extension of the low-energy bremsstrahlung components that we deter- 
mined from fits to the GBM data; it must be from a distinct population of electrons extending 
to energies of several hundred MeV. We have determined the power-law parameters (PL3) 
of the fitted bremsstrahlung spectrum and list them in Table 1. However, this third high 
energy electron component cannot be produced by the acceleration models mentioned §4.1 
above, which produce spectra that steepen beyond t ens of MeV due to synchrotron energy 
losses that increase with energy (see lPark et al.lll997l ). and must have a quite different origin. 
Consequently we believe that this is a less likely scenario than the hadronic model. 

Assuming that the LAT emission is from hadronic interactions, we have fit the LAT 
spectrum with calculated pion-decay spectra produced by accelerated ions having differential 
power-law indices from —2.5 to —7.5. With 67% confidence (based on x^) we conclude that 
the spectrum of accelerated ions responsible for the pion-decay emission must be steeper 
than a power-law with index —4.5. We note that there is no change in the quality of the 
fits for indices steeper than —5 due to limited statistics >400 MeV. We list the fluence of 
pion-decay photons >200 keV derived from these fits in Table 1. This fluence is a factor of 
ten below the limit derived from fitting the GBM BGO data up to 40 MeV, illustrating the 
significantly greater sensitivity of the LAT instrument. 

We can use the results of our GBM and LAT spectral anal yses to obtain information 
on ions accelerated in the impulsive phase of the June 12 flare. iMurphy et al.l (119971 ) have 
described how parameters derived from integrated spectroscopic fits and temporal studies 
can be used to obtain this information. We first use the nuclear de-excitation line, neutron- 
capture line, and pion-decay fiuences listed in Table 1 to estimate the overall shape of the 
accelerated ion spectrum. These three emissions are produced by accelerated ions within 
distinct energy ranges: ~5-20 MeV for the de-excitation lines, ~ 10-50 MeV for the neutron 
capture line, and >300 MeV for the pion-decay emission. Ratios of these emissions therefore 
determine the relative numbers of accelerated ions in the associated energy ranges. We then 
obtain spectral indices across t hese energy ranges by cornparing measured ratios with ratios 
from theoretical calculations (IMurphy et al.l 119871 . l2005l . 120071 ) based on updated nuclear 
cross sections. 



If we assume that the LAT emission >30 MeV was entirely due to pion-decay emission. 
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then we estimate that the flare-accelerated ion spectrum was consistent with a series of power 
laws, softening with energy, with indices of ~— 3.2 between ~ 5 — 50 MeV, ~— 4.3 between 
~50-300 MeV, and softer than ~— 4.5 above 300 MeV. However, these calculations assume 
that these individual power laws continued without break to high energies, which cannot be 
the case where the spectrum softens with increasing energy. We will describe a more refined 
representation of the accelerated proton spectrum in a forthcoming publication. 



4.3. Combined GBM and LAT Timing Studies 

The combined GBM and LAT observations also provide us with the opportunity to 
study acceleration and transport in this impulsive flare. In Figure [5] we plot 5-s resolution 
time profiles of the fitted bremsstrahlung and nuclear de-excitation line fluxes from the GBM 
and the >100 MeV flux observed by the LAT. For purposes of comparison we have plotted 
the bremsstrahlung profile over the nuclear and LAT >100 MeV histories in b) and c), 
respectively. We have also plotted the nuclear profile over the LAT >100 MeV history. The 
early peaking of the bremsstrahlung suggests that the higher-energy emissions were delayed 
by a few seconds and a cross-correlation study indicates that the overall lag between the 
LAT >100 MeV flux and the GBM 300 keV bremsstrahlung flux is ~3 s. 

This delay warranted a higher-time resolution comparison of 100 - 500 keV GBM and 
LAT >30 MeV counting rates rather than fitted fluxes. The 100 - 500 keV band is dominated 
by electron bremsstrahlung as can be seen in Figure [31 In Figure [6] we plot GBM/BGO 
100-500 keV rates at 320 ms resolution with 3-s LAT/LLE >30 MeV rates over plotted. 
The hard X-ray profile reveals the presence of a clearly separated initial peak along with 
other structures. The onset of the > 30 MeV emission appears to be ~3 s following the 
bremsstrahlung and rises to a peak about 10 s after the 100 -500 keV peak. The LLE 
profile appears to reflect the double-peaked bremsstrahlung profile with a delay of about 10 
seconds. From a cross correlation analysis of the two profiles plotted in Figure [6] we find 
that the >30 MeV emission lags the bremsstrahlung by 6 ± 3 s. There are two fundamental 
implications of the time profiles in Figure |6l 1) protons and/or electrons began reaching 
energies above 100 MeV within a few seconds of the time it took to accelerate electrons to 
energies of hundreds of keV; and 2) the overall acceleration time scale of the >100 MeV 
particles is similar to that observed in hundreds of keV electrons, but delayed by about 10 
seconds. 
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Table 1. Best- fitting spectral parameters 



Pcirameter Value 

PLl fluence at 300 keV 2.85 ± 0.1 7 cm-^ keV-^ 

PLl index 3.31 ± 0.09 

PL2 fluence at 300 keV 0.08 ± 0.02 7 cm^^ keV"! 

PL2 index < 1.2 

PL2 Exponential Energy 2400 it 800 keV 

0.511 MeV line fluence'' 11.3 ± 2.5 7 cm-^ 

2.223 MeV line fluence* 21.3 ± 2.0 7 cm-^ 

Nuclear line fluence 23.5 ± 2.5 7 cm~^ 

Pion-dccay fluence (GBM) > 200 keV 1.5 ± 2.5 7 cm-^ 

Pion-decay fluence (LATt- ) > 200 keV 0.62 ± 0.07 7 cm-2 

Pion-decay fluence (LAT) > 100 MeV 0.13 ± 0.015 7 cm-^ 

PL3 fluence at 30 MeV (9.2 ± 2.0) xlO-'^7 cm-^ keV"! 

PL3 index 1.9 ± 0.2 

'^Integrated from 00:55:40 to 00:59:50 UT 

''Computed by extrapolating to low energies the model that best flts LAT 
data 
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Fig. 5. — Time profiles of the a) bremsstrahlung differential flux at 300 keV (GBM), b) total 
nuclear de-excitation line flux (GBM), and c) integral pion-dccay flux >100 McV (LAT). 
Solid curves over-plotted in b) and c) are the arbitrarily normalized bremsstrahlung flux 
plotted in a) ; dashed curve over-plotted in c) is arbitrarily normalized nuclear line flux from 
panel b). 
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4.4. Search for >100 MeV Emission Following the Flare 



The Fermi LAT detects quiescent emission from the Sun on a near daily basis ( lAbdo et al 



201ll ). This emission comes from cosmic-ray proton interactions in the solar atmosphere 
and photosphere, and from Compton scattering of cosmic-ray electrons on solar black- 
body photons. The LAT is therefore a sensitive monitor of temporally extended solar-flare 
emission such as detected by CGRO E GRET experiment following the 1991 June 11 flare 
(IKanbach et al.l Il993l : iRank et al.ll200ll ). We therefore studied the emission within 15° of 
the Sun in the hours preceding and following the flare. The standard LAT data products 
were used in the analysis, which modele d the region around the Sun includin g all sources in 
the Second Fermi LAT Source Catalog (IThe Fermi-LAT Collaboratiorul201ll ). isotropic and 
Galactic diffuse emissions, and spatially extended Compton-scattered solar photons discussed 
above. 

In Figure [7] we plot 95% confldence limits on the >100 MeV flux from the solar disk in 
~30 min exposures every two orbits from 6 hr before the flare to 22 hr after it. We note 
that the Sun was outside the FOV for LAT standard-product analysis during the flare, but 
its 7-ray emission could be studied in the orbits just before and after the flare. During this 
time period the Moon passed within 10° of the Sun and its flux of ~ 1 x 10^^ 7 cm^^ s^^ 
made a signiflcant contribution to the measured solar fluxes because we did not include it 
in the model of the source region. There is no evidence for an increase in the solar emission 
following the time of the flare, which is denoted by the peak flare flux >100 MeV. This peak 
flux is about 1000 times higher than the plotted upper limits. We wish to compare this to 
the time extended >100 MeV emission observed by EGRET following the 1991 June 11 flare 
(IKanbach et al.lll993l : iRank et al.ll200ll ). But EGRET was saturated at the peak of the flare 
and therefore its time history cannot be normalized to the >100 MeV LAT peak flux of the 
June 12 flare. However, after the peak of the flare the ratio of the >100 MeV EGRET flux 
to the 2.223 MeV neutron-capture line measured by COMPTEL, which operated normally 
throughout the June 11 flare, was relatively constant. We therefore used the COMPTEL 
time history as a proxy for the EGRET >100 M eV photon fluxes . The time proflle of the 
COMPTEL 2.223 MeV neutron-capture line flux (jShare et al.lll997l ). normalized to the peak 
>100 MeV LAT June 12 flux, is shown in the Figure [71 We see that during the first 30 
minute exposure following the flare, the LAT upper limit is a factor of ~ 20 below the value 
expected if there had been comparable extended emission similar to that found on 1991 June 
11. 
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Fig. 6. — Comparison of high-time resolution profiles of the 100 - 500 keV emission ob- 
served in the GBM BGO detector and of >30 MeV LLE data. A cross- correlation analysis 
indicates that the high-energy 7-ray emission had an overall lag of 6 ± 3 s relative to the 
bremsstrahlung. 
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Fig. 7. — 95% confidence upper limits of >100 MeV solar 7-rays measured by LAT within 7 
hours preceding and 22 hours following the 2010 June 12 flare. Dashed line is the 1991 June 
11 time history of the 2.223 MeV neutron-capture line observed by COMPTEL experiment 
and diamonds are EGRET >50 MeV data both normalized to the peak >100 MeV flux 
measured by the LAT for the June 12 flare. The 2010 June 12 flare does not show extended 
emission like the 1991 June 11 flare observed with EGRET. 
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5. Summary and Discussion 

The 2010 June 12 flare was the flrst in Cycle 24 to be observed to emit nuclear 7 
rays. It was also the flrst flare detected by the Fermi LAT at energies above 30 MeV. The 
hard X-ray and nuclear line radiation was observed both by the Fermi GBM and RHESSI 
spectrometers. In this paper we only analyzed GBM data because RHESSI was offset from 
the Sun to study the Crab Nebula during the time period of the flare; this affected knowledge 
of the instrument response. 

The fact that the flare emitted detectable 7 rays at al l is surprisi n g bec ause its peak soft 



X-ray emission only reached a GOES M2 level. However, IShih et al.l ( l2009l ) and others have 
shown that 7-ray line fluences are only weakly correlated with GOES soft X-ray emission 
but are strongly correlated with electron bremsstrahlung fluences >300 keV. This is true for 
the June 12 flare as we flnd that the measured bremsstrahlung and 2.223 MeV fluences are 
consistent with the established correlation. 

The flare originated from a compact region and its hard X-ray emission only lasted 50 
s. Figure 6 reveals striking information about the processes that accelerate protons and/or 
electrons to energies of hundreds of MeV. We flnd that although some of the particles reach 
energies ^100 MeV within about 3 s, the bulk of these particles reach such high energies 
following a delay of about 10 s. This is revealed in the delayed double-peaked time structure 
>30 MeV that is similar to what is observed in hard X-rays. 

In Table 1 we list the best-fltting parameters from our flts to the GBM and LAT spectra. 
These include the amplitudes (at 300 keV) and indices of two power-law continua observed by 
the GBM between 300 keV and 8 MeV. The flrst is an extension of the hard X-ray spectrum 
observed by the GBM Nal detectors. The second appears to be a hard power law with an 
exponential cutoff energy near 2.5 MeV. Although the GBM only has moderate spectral 
resolution, it was able to measure the fluences of the 0.511 MeV annihilation and 2.223 MeV 
neutron capture lines, and the total nuclear de-excitation emission. There was no evidence in 
the GBM data for flare emission above about 8 MeV. In contrast the LAT detected signiflcant 
continuum emission from ~30 to 400 MeV about an order of magnitude below the GBM 
upper limits. This reflects the excellent sensitivity of the LAT for observing solar flares. 
This radiation could be either pion-decay emission or primary electron bremsstrahlung. In 
Table 1 we list the measured pion-decay fluences >200 keV and >100 MeV. We also list the 
parameters of a high-energy power-law bremsstrahlung component that flts the LAT data 
equally well. Theoretical arguments made in §4.2 weigh against an electron origin for the 
emission observed by the LAT. Under the assumption of an hadronic origin for the LAT 
emission and using 7-ray line measurements by the GBM, we have estimated the shape of 
the accelerated-ion spectrum that could have produced the combined spectrum. 
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We have also set significant constraints on any time-extended >100 MeV emission. Our 
limit on the >100 MeV flux of photons is about an order of magnitude below what would 
have been expected if the decay followed th at observed in the well-studied 1991 June 11 
flare (IKanbach et al.lll993t iRank et al .1120011 ). Thus, there is no evidence for precipitation 
of trapped flare particles, particles accelerated in magnetic loops after the impulsive phase, 
particles accelerated in CME-associated reconnection sheets (lRyanll2000l ). or p articles sharing 



the same origin as the Solar Energetic Particles (SEPs) observed in space (IRamaty et al. 



1987 



Oliver et al. 



1993h. We estimate from white light measurements of the 1991 June 11 
flare ( ISakurai et al.l Il992l ) that the foot-point separation was ~ 2.5 x 10^ km, about 2.5 



times larger than the June 12 flare. It is possible that longer coronal loops are necessary for 
time-extended acceleration and/or trapping of protons. 

This paper primarily addresses only the 7-ray observations of the 2010 June 12 flare. 
Other studies are currently in progress involving hard X-ray observations of the flare by the 
Fermi GBM Nal detectors and comparisons of the characteristics of the ion and electron 
populations at the Sun and in space and at Earth. 
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